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Abstract 

We find that there is no supersymmetric flavor/CP problem, /x-problem, cosmological 
moduli/gravitino problem or dimension four/five proton decay problem in a class of su- 
persymmetric theories with 0(1) GeV gravitino mass. The cosmic abundance of the non- 
thermally produced gravitinos naturally explains the dark matter component of the universe. 
A mild hierarchy between the mass scale of supersymmetric particles and electroweak scale 
is predicted, consistent with the null result of a search for the Higgs boson at the LEP-II 
experiments. A relation to the strong CP problem is addressed. We propose a parametriza- 
tion of the model for the purpose of collider studies. The scalar tau lepton is the next to 
lightest supersymmetric particle in a theoretically favored region of the parameter space. 
The lifetime of the scalar tau is of 0(1000) seconds with which it is regarded as a charged 
stable particle in collider experiments. We discuss characteristic signatures and a strategy 
for confirmation of this class of theories at the LHC experiments. 



1 Introduction 



In spontaneously broken supersymmetric theories, there is a spin-half Goldstino fermion which 
is eaten by the gravitino as its longitudinal components. By supersymmetry, the Goldstino 
must be accompanied with its superpartner whose spin is zero if supersymmetry is broken by a 
vacuum expectation value of the F-component of a chiral superfield. A chiral supermultiplet is 
formed by the Goldstino, its scalar superpartner, and the non-vanishing F-term, which we call 
the chiral superfield S. The low energy physics is then described by matter superfields, gauge 
superfields and the chiral superfield S. 

There are variety of possibilities for couplings between matter /gauge superfields in the 
supersymmetric standard model and the superfield S. These possibilities have been classified as 
follows. If we assume that the couplings are suppressed by the Planck scale (Mpi), such as £ 9 
[{S/Mpi)W'^Wa]e2 with TV" being gauge fields, the model is called the "gravity mediation" |]l|, ^. 
Another possibility that the gauge kinetic function is of the form, C 3 [(logS'/(47r)^)VF"Wa]02, 
is called the "gauge mediation" 1^, ^, ^, ^. This is the form we obtain after integrating out 
vector- like fields which obtain masses proportional to (5) [0]. If the coupling is more suppressed 
than the Planck scale, effects of the "anomaly mediation |^" give the largest contribution to 
supersymmetry breaking terms in the Lagrangian. Among those scenarios, gauge mediation 
assumes the strongest interaction between the matter /gauge fields and S while the anomaly 
mediation effects are the weakest. The size of supersymmetry breaking, Fs, therefore has a 
relation, -Pgauge -^gravity ^ -^anomaly) when we fix the scale of gaugino/sfermion masses. The 

.,• gauge gravity anomaly 7^ 

gravitmo masses are rriy^ ^ "^3/2 "^3/2 ^ ^3/2 ^ ^S- 

The question is what size of the gravitino mass (i.e., the supersymmetry breaking scale) is 
preferred by phenomenological and cosmological requirements. This is an interesting question 
since each scenario predicts a quite different pattern in the spectrum of the supersymmetric 
particles, which we will search for at the LHC experiments. Strategies for finding supersymmetric 
particles and measurements of model parameters will be also different for different scales of 7713/2- 

There have been many model-building efforts in making supersymmetric models realistic 
in each category: gravity, gauge or anomaly mediation. If one of them had been completely 
successful, we could have believed in the scenario and used the model as the standard super- 
symmetric model. However, unfortunately, there is no such a standard model so far because 
of the fact that neither of these scenarios are fully realistic by different reasons. In the gauge 
and anomaly mediation scenarios, there is a problem with the electroweak symmetry breaking, 
i.e., the //-problem. The pure anomaly mediation, in addition, predicts tachyonic scalar leptons 
which are not acceptable. Although the gravity mediation scenario does not suffer from those 
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problems, it has been known that sizes of flavor and CP violation are expected to be too large. 
There are also cosmological constraints. In particular, in gravity mediation models, a moduli 
problem caused by fields in the supersymmetry breaking sector destroys cosmological successes 
of the (supersymmetric) standard model such as the big-bang nucleosynthesis (BBN) and 
also cold dark matter by thermal-relic neutralinos. 

In this paper, we reconsider problems in supersymmetric models by using an effective field 
theory described by the field S and the matter /gauge fields. By doing so, we can discuss each 
of these scenarios as a different choice of functions of S which define an effective theory. The 
labeling can be done by projecting the function space onto a one-dimensional axis of the gravitino 
mass. In this formulation, we find that there is a sweet spot in between the gauge and gravity 
mediation ~ GeV) where the theory is perfectly consistent with various requirements. 

All the classic problems, such as the flavor/CP problem and the /i-problem are absent. The 
theory also avoids a cosmological moduli problem caused by the scalar component of S. Non- 
thermally produced gravitinos through the decay of the 5-condensation naturally account for 
dark matter of the universe. A simple ultraviolet (UV) completion of the theory exists, which 
is actually a model of grand unification without neither the doublet-triplet splitting problem 
nor the proton decay problem. Relations to the strong CP problem and the supersymmetric 
fine-tuning problem are also addressed. We discuss a characteristic spectrum of supersymmetric 
particles, and demonstrate how we can confirm this scenario. 

In the next section, we rewrite the various supersymmetric models in terms of the effective 
Lagrangian described by the Goldstino multiplet S and particles in the minimal supersymmetric 
standard model (MSSM). The section includes review of the supersymmetry breaking and its 



transmission. A concrete set-up is defined in subsection 2.4 and discuss its successes there. We 
then discuss low energy predictions of the framework in Section ^ A parametrization of the 
model and a way of calculating the spectrum of supersymmetric particles are presented. Collider 
signatures are discussed in Section ^ We demonstrate a method of extracting model parameters 
in the case where the scalar tau (stau) is the next to lightest supersymmetric particle (NLSP). 



2 Theoretical set-up 



We construct a phenomenological Lagrangian of the supersymmetric standard model and con- 
sider various requirements from particle physics and cosmology. We will arrive at a scenario 
with ~ 1 GeV. 
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2.1 S sector 



We derive here a description of a super symmetry breaking sector by the Goldstino chiral 
superfield S. This corresponds to the construction of the Higgs sector in the standard model. 
As any models of the electroweak symmetry breaking flow into the standard model with various 
mass ranges of the Higgs boson at low energy, the model below provides a standard low energy 
description of a variety of supersymmetry breaking models. 

We concentrate on F-term supersymmetry breaking scenarios as most of the supersymmetry 
breaking models are of this type. To ensure a non-vanishing vacuum expectation value of the 
F-component of a chiral superfield S, we add a source term in the Lagrangian: 

jC 3 n?Fs + h.c. (1) 

This term can be expressed in terms of the superfield as follows: 

W B m^S . (2) 



We can also write down an arbitrary Kahler potential, Kg, for the kinetic and interaction 
terms of S. As long as d'^Ks/{dSdS^) is a non-singular function, Fs ^ is obtained by the 
equation of motion. For example, the low energy effective theory of the O'Raifeartaigh model [0] 
has a Kahler potential: 

Ks = S^S-^-^, (3) 

where A is the mass scale of the massive fields which have been integrated out. In general, if 
S carries some approximately conserving charge, the Kahler potential is restricted to the form 
in Eq. (|^) (up to the sign of the second term).|^ The second term gives a mass to the scalar 
component of 5, ms- 

,2 



2Fs 2m2 /- /MpA 
ms = — = ^ = 2V3my, ^-^ ) 



(4) 



and stabilizes the value of S at 



5 = 0. (5) 



Here we ignored supergravity effects. The fermionic component of S remains massless. This is 
the Goldstino fermion associated with the spontaneous supersymmetry breaking. 



*In fact, the cubic term in the Kahler potential, K B +h.c. can be eliminated by the shift of S in general. 
However, once we take into account interaction terms between S and the MSSM fields, the origin of S has a 
definite meaning and we cannot shift away the cubic term. 
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Note that the existence of the chiral superfield S in the above effective theory does not nec- 
essarily mean that the supersymmctry breaking sector contains a gauge singlet chiral superfield 
in the UV theory. The S field can originate from a component of some multiplets or can be a 
composite operator in physics above a 'cut-off' scale A. It is totally a general argument that there 
is a gauge singlet chiral superfield S in the effective theory below the scale of supersymmetry 
breaking dynamics, A, as long as A^ > m?. 

The Lagrangian discussed above is analogous to the Higgs sector in the standard model. The 
two parameters and m^/A^(~ nig/m?) correspond to the parameters and A//(~ mf^/v'^) 
in the Higgs potential, V = (A///4)(|i/p — f^)^. We should not trust this effective theory if 

I !^ ^ \/47r as it violates the unitarity of scattering amplitudes of the gravitinos at high 
energy just like the standard model with Xh > An. 

2.2 Matter /gauge sector 

The superpotential of the MSSM is 

H^MSSM = QHuU + QHdD + LHdE + fiH^Ha , (6) 

where we suppressed the Yukawa coupling constants and flavor indices. The last term, the fi- 
term, is needed to give a mass to the Higgsino, but it should not be too large. For supersymmetry 
to be a solution to the hierarchy problem, i.e., ^ Mp\, the /i-term is necessary to be of 

the order of the electroweak scale (or scale of the soft supersymmetry breaking terms) . 

This is called the ^-problem. The fact that ^ is much smaller than the Planck scale suggests 
that the combination of HuH^ carries some approximately conserving charge. 

There are many gauge invariant operators we can write down in addition to the above 
superpotential such as 

Wiji = UDD + LLE + QLD , (7) 

and 

Wdim.5 = QQQL + UDUE . (8) 
These are unwanted operators as they cause too rapid proton decays. 

The yU-problem and the proton decay problem above are actually related, and there is a 
simple solution to both problems. The Peccei-Quinn (PQ) symmetry with the following charge 
assignment avoids too large /x-term and the proton decay operators. 

PQ(Q) = PQ{U) = PQ{D) = PQ{L) = PQ{E) = - ^ , (9) 
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PQ{H^) = PQ{Hd) = 1 



(10) 



This symmetry is broken explicitly by the //-term, PQ{p) = —2. Since it is a small breaking 
of the PQ symmetry, the coefficients of the dimension five operators are sufficiently suppressed. 
The unbroken Z4 symmetry, which includes the i?-parity as a subgroup, still forbids the su- 
perpotential terms in Eq. (0) and ensures the stability of the lightest super symmetric particle 
(LSP), leaving us to have a candidate for dark matter of the universe. 

The Majorana neutrino mass terms, W B LLHuHu, are forbidden by the PQ symmetry, 
but large enough coefficients can be obtained by introducing another explicit breaking of the 
PQ symmetry. For example, we can write down LLHuHu/Mj^j with PQ{M]\f) = 1 without 
introducing proton decay operators or too large /x-term. The Z4 symmetry above is broken 
down to the i?-parity with this term. 

In fact, there is another symmetry which can play the same role as the PQ symmetry, called 
i?-symmetry. The charge assignment is 

R{Q) = R{U) = R{D) = R{L) = R{E) = 1 , (11) 

R{Hu) = R{Hd) = . (12) 

Again, R{ii) = 2 explicitly breaks the i?-symmetry down to the i?-parity. In this case, the 
LLHuHu term is allowed by the symmetry. 

In summary, there are approximate symmetries, U(1)pq and U(l)/j, in the Lagrangian of 
the MSSM. If one of them is a good (approximate) symmetry of the whole system, it provides 
us with a solution to the /z and the proton decay problems. 

2.3 Interaction to mediate the supersymmetry breaking 

Now we discuss interaction terms between the S-sector and the MSSM sector. These interactions 
determine the pattern of supersymmetry breaking parameters which are relevant for low energy 
physics. We review here three famous mechanisms; gravity, gauge, and anomaly mediation 
models, as choices of the form of the interactions. Each of these scenarios suffer from different 
problems. Understanding nature of those problems guides us to a phenomenologically consistent 
model. 
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2.3.1 Gravity mediation 



The simplest scenario is to assume general interaction terms suppressed by the Planck scale. 
This is called the gravity mediation. The Kahler potential is 

^gravity- ^2^ + ^ ^pi + ^•"'y' + ^^"^^ 



s^s{hIh^ + hIh,) 

where $ represents the quark and lepton superfields in the MSSM, and we omit 0(1) coefficients. 
Planck suppressed operators in the gauge kinetic function generate gaugino masses: 



/gravity =(^ + 4^^) W-W^ , (15) 



where g is the gauge coupling constant. 

The first term in Eq. ( [T^ ) and the second term in the bracket in Eq. (^5|) generate sfermion 
masses and gaugino masses, respectively. Both of them are of 0{Fs/Mpi) ~ Therefore, 
the gravitino mass, m^/2j is O(IOO) GeV in this scenario. The /i-problem is completely solved 
in a quite natural way [0]. The first and second terms in Eq. (|l^ ) generates /u ~ O (7713/2 )• 

This mechanism for the /x-term generation is consistent with the discussion in the previous 
subsection. The i?-symmetry introduced before can be preserved once we assign R(S) = 0. 
The term in Eq. (|2|) breaks the /^-symmetry by R{m?) = 2 at the intermediate scale, m? = 
\/37773/2Afpi. This is still small enough for the proton decay operators. 

On the other hand, the natural solution to the /x-problem is not compatible with the PQ 



symmetry. The term in Eq. (15) restricts the PQ charge of S to be vanishing, whereas the 
term responsible for the //-term generation, K 3 HuH^, determines that PQ{S) = 2. The 
PQ symmetry must, therefore, be maximally violated. This implies that none of the terms 



in Eq. (14) can be forbidden by approximate symmetries of the theory. This fact becomes 
important in the discussion of the super symmetric CP problem. 

Even though the //-problem is solved perfectly, there are several serious problems in this 
scenario. Since there is no reason for the alignment of the flavor structure in the first term in 
Eq. (|T^), too large rates for flavor changing processes are predicted. We expect flavor mixings 



of 0(1) from this form of Lagrangian. Such large mixings are unacceptable unless the sfermion 
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masses are of 0(10) TeV or heavier [|T2|. The CP violating phases in the supersymmetry breaking 
terms are also expected to be 0{1). In particular, a phase of the combination, mi/2fJ'iB^)* , 
with mi/2 the gaugino mass and Bfi defined by £ 9 B^HuHd + h.c, cannot be eliminated 
by field redefinitions. The /U and terms are generated from the multiple terms in Eq. (|l^ ) 
with different weights, leading to non-aligned phases generically. With an 0(1) phase for the 
combination, constraints from the electric dipole moment of electron, for example, push the 



mass limits of supersymmetric particles to be 0(10) TeV [12|. 

There is another serious problem in cosmology. Due to the terms in Eq. (p^), the scalar 
component of S cannot carry any (even approximately) conserving charge. In this case, there 



is a moduli problem y, 13]. The value of S after the inflation is displaced from the minimum 
due to the deformation of the S potential during inflation, and at a later time S finds its true 
minimum and starts coherent oscillation about the true minimum. The energy density of the 
oscillation then dominates over the universe unless the displacement is much smaller than the 
Planck scale. The decay of S, in turn, either destroys the success of the BBN @ or overproduce 
gravitinos [^] depending on the mass range of ms (see [14, ^] for earlier works). There is no 
range of ms which is consistent with the cosmology |T^ ]. 

It has been widely accepted that the lightest neutralino accounts for dark matter of the 
universe in gravity mediation scenarios. Abundance of thermally produced neutralinos can be 
calculated without information on the detail history of the universe, and we obtain the correct 
order of magnitude. However, once we take into account the existence of the superpartner of 
the Goldstino, S, which always exists, the successful cosmology is spoiled. We argue that an 
assumption made in the standard calculation that the universe was normal up to temperatures 
of O(IOO) GeV is inconsistent with the structure of underlying models. 

2.3.2 Gauge mediation 

Some of shortcomings in gravity mediation can be cured in gauge mediation models. We assume 
in gauge mediation that the size of the supersymmetry breaking, Fs (and therefore ?ti-3/2)5 is 
much smaller than that in gravity mediation. The contributions to the soft supersymmetry 
breaking terms come from 

^fcr^ = -^^yF^2(i?)(log |5|)2$t$ , (16) 
1/1 2iV--s i^g5W«^^ ^ (17) 
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where C2{R) is the quadratic Casimir factors for fields ^. These terms are generated by 
integrating out A'mess numbers of messenger fields, / and / in the fundamental representation, 
which have couplings to S in the superpotential, W 3 kSff 0] . Singularities at S = indicate 
that the messenger fields become massless at the point .|^ Therefore, the theory makes sense only 
if the potential of S has a (local) minimum at 7^ 0. Low energy parameters depend on the 
coupling constant k only through a logarithmic function. The dependence is encoded as the 
messenger scale Mmess = k{S) at which gauge mediation effects appear. 

The contributions from gauge mediation are much larger than those from gravity mediation 
in Eqs. (1^,14,15) provided that the value of S is stabilized at 5 <C Mpi. Since there is no flavor 



dependent terms in Eq. (|16|), due to the flavor blindness of the gauge interactions, constraints 
from flavor violating processes can be easily satisfied when ^ 0(1) GeV. 

Another interesting feature is the enhancement of the S couplings to the MSSM particles 
The scalar component of S now has couplings to gauginos, A: 

£9^5AA + h.c., (18) 

which can be much larger than the coupling to gravitinos, V'3/2- 

C 3 -^SH^/2i^^/2 + h.c. , (19) 



depending on the value of {S). Therefore, the branching fraction of the S decay into gravitinos is 
suppressed and the situation of gravitino overproduction from the S decay can be ameliorated. 

However, unfortunately, by lowering the gravitino mass 7713/2, we have lost the natural 



mechanism for generating a //-term. The contributions from Eq. (14) to the /i-term are too 
small. It is possible to obtain a correct size of //-term by assuming a direct coupling between S 
and Higgs fields such as 

W 3 eSHuHd, (20) 

with a small coefficient e |p. This term, however, predicts too large term, ~ (47r)^m]^/2i 

which is unacceptable from the electroweak symmetry breaking. The conclusion is the same if 
we try to generate a /x-term from Kahler terms, e.g., 

K^J^2^H^H, + \..c. (21) 



^In this discussion, we have defined the origin of S to be the point where the messenger particles become 



massless. It is not necessarily the same definition in subsection 2.1, We will discuss a whole set-up together with 
the S'-sector shortly. 
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This term can be generated by integrating out messenger fields if there is an interaction between 
the Higgs and the messenger fields [18, 19 j. Although this term induces the correct size of 
/i-term, /i ~ mi/2) Bfi/ fx is again larger than mi/2 by a one-loop factor]^ 

The /i-problem in gauge mediation models cannot be solved by going to the next to minimal 
super symmetric standard model (NMSSM). The correct electroweak symmetry breaking is not 
achieved without further extensions of the model pO| ]. 

We cannot discuss the supersymmetric CP problem without specifying the mechanism for 
//-term generation because the physical phase arg(?7ii/2/^(^/^)*) is not determined. 

Although it is slightly model dependent, there is another issue in gauge mediation models. 
In many supersymmetry breaking models, S carries a conserving charge. For example, in the 
O'Raifeartaigh model, there is an unbroken i?-symmetry where S carries charge 2. In this case, 



as discussed in subsection |2.1| , the S field is stabilized at 5 = where we cannot integrate out 
the messenger fields (it is the singular point of the effective Lagrangian) Additional model 
building efforts to shift the minimum of the S potential have been needed in this type of models. 

More explicitly, what we need is to spontaneously or explicitly break the i?-symmetry in 
supersymmetry breaking models. If we break it explicitly in the Lagrangian, the theorem of [^] 
says that a supersymmetric minimum appears somewhere in the field space. Recently, there 
have been extensive studies on this subject, and many simple models with explicit breaking 



of ii-symmetry have been proposed |22, 23, 25 1 by allowing a meta-stable supersymmetry 



breaking vacuum |26|. (See also |2^, ^] for recent models with spontaneous breaking of 
i?-symmetry.) An obvious possibility is to add an i?-breaking cubic term in Eq. (^) with a small 
coefficient e, 

5K = ^S^S"^ + h.c. (22) 
A 

This term shifts the minimum of S* to 5 = eA/2. This is equivalent to give a small mass term to 
the messenger fields W B eAff [^] by a field redefinition {S — eA/2) — > 5. A small mass term 
for S, W B eS^ also shifts the minimum of S. 

In fact, it has been known that these ad hoc deformations of the model were not necessary 



* There is a logical possibility of generating terms like 

^ J^H^Hd\og\S\+h.c. , 

if the Higgs fields have interactions with messenger fields. In this case, the B/i-term is not generated, whereas the 
/i-term is generated with the same size as the gaugino masses. This is perfectly consistent with the electroweak 
symmetry breaking and also the absence of the CP phase in mi/2/i(-B/i)* as Bjj, = 0. This is the only term which 
can be written down if S carries an approximately conserving charge. However, the authors are not aware of an 
explicit model to realize this situation. 

^The origin of S is now uniquely determined once we assign a charge to S. 
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once we take into account supergravity effects |22|. The gravity mediation effects generate a 
linear term of S in the potential, 

V B 2m3/2m25 + h.c. (23) 
This is a soft supersymmetry breaking term associated with the linear term in the superpotential 

41 



in Eq. (^). By balancing with the mass term, V 3 m'g\S\'^ with ms in Eq. (^), we obtain 



This shift is due to the fact that i?-symmetry must be broken explicitly in the supergravity 
Lagrangian (by the constant term in the superpotential) in order to cancel the cosmological 
constant ||l^. By taking a large A, the shift can be arbitrarily large. Note here that the shift 
is not suppressed by the gravitino mass which characterizes the effects of gravity mediation. 



This phenomenon has been known as the tadpole problem for singlet fields ||2|, 30, 31 1. Small 
soft supersymmetry breaking terms destabilize the hierarchy if there is a gauge singlet field. 
However, this is not a problem at all for the field S and it is even better in gauge mediation 
to have a large enough vacuum expectation value of S. Since this effect always exists, it is the 
most economical way of having S ^ 0. 



2.3.3 Anomaly mediation 

If there is no direct coupling between S and the MSSM particles even including Planck scale 
suppressed operators, the leading contribution to the sfermion/gaugino masses comes from 
anomaly mediation effects: 

a'^b 2 1 2 /OtrN 

^i/2 = j^my,, = -^^^m3/2 , (25) 

where b and 7 are the beta function coefficient and the anomalous dimension, respectively, and 
is the renormalization scale For having = O(IOO) GeV, a large gravitino mass 

~ 10 — 100 TeV is needed. There are several good features of this scenario. Because 
of flavor blindness of the mediation mechanism, there is no super symmetric flavor problem. 
The large value of enhances the decay rate of the gravitino. This makes the gravitino 

cosmologically harmless as it decays before the BBN starts. The cosmological moduli problem 
is also absent. The S field can have any conserving charges, and thus it is reasonable to assume 
that S has stayed at the symmetry enhanced point, S = 0, during and after inflation so that 
there is no large initial amplitude. 

Unfortunately, the minimal model is inconsistent with the observation. The scalar leptons 
have tachyonic masses which would cause a spontaneous breaking of U(l)em and makes the 
photon massive. Therefore, we need a modiflcation of the model. 
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Also there is a yU-problem which is very similar to the situation in gauge mediation. One 
can assume small couplings between S and the Higgs fields to give a /i-term, but it causes a too 
large Bfi-teim, -B/x/// ~ ^-3/2 with 7713/2 ~ 10 — 100 TeV which is unacceptable. 

2.4 Sweet Spot Supersymmetry 

We have encountered many problems in gravity, gauge and anomaly mediation models. Those 
are summarized as follows: 



• Gauge mediation (777,3/2 ^ 100 GeV) 
Problems: fi, (CP), 

• Gravity mediation (7773/2 ~ 100 GeV) 
Problems: Flavor, CP, moduli, 

• Anomaly mediation (7773/2 ~ 10 — 100 TeV) 
Problems: fi, tachyonic sleptons, (CP). 

There have been many attempts to circumvent these problems. For example, in Ref. ^] 
it has been proposed to extend a model of gauge mediation to the NMSSM by introducing a 
new singlet field. However, for the successful electroweak symmetry breaking, further extension 
of the model were necessary such as introduction of vector-like matters. Similar attempts have 



been done in Ref. 1^^, 33] in anomaly mediation models. The gaugino mediation [^] is a variance 
of the gravity mediation and known to be a successful framework for solving the flavor problem. 
However, since the model relies on the SW^Wa term for the gaugino masses, the moduli problem 
and the CP problem remain unsolved. In Ref. |^5|, a mixture of anomaly and gauge mediation 
is proposed as a solution to the tachyonic slepton problem (see also [^]). The idea is to modify 
the structure of the anomaly mediation by introducing an additional light degree of freedom, X. 
It is claimed that the /^-problem and the tachyonic slepton problem can be solved by assuming 
appropriate couplings of X to the messenger and Higgs fields. However, it is unclear whether 
such a light degree of freedom is consistent with cosmological history. 

It is interesting to notice here that gauge and gravity mediation scenarios do not share 
problems. This fact motivates us to think of theories in between gauge and gravity mediation. 
The idea is to solve flavor and moduli problem by reducing 7773/2, and solve the /i-problem in a 
similar fashion to the gravity mediation models. The CP problem can also be solved because we 
can have an approximate PQ symmetry to forbid the Bfi-term so that arg(777^/2/i(i?/u)*) = 0. 
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Figure 1: Schematic picture of mediation mechanisms. Different mechanism works for different 
values of gravitino masses. A sweet spot exists at m3/2 ~ 1 GeV where there is no 
phenomenological or cosmological problem. 

The sweet spot exists at ra-^i2 ~ 1 GeV (see Fig. [l|). The interaction terms between 
matter /gauge field and S are the same as those in gauge mediation (Eqs. ( p^ , p!7[ ) ) . For ra^^j^ ~ 
1 GeV, possible flavor violating contributions from gravity mediation in Eq. ( p!3| ) are sufficiently 
small. The couplings of S to Higgs fields are 

Here we replaced the Planck scale in Eq. (|l^ with the "cut-off" scale A introduced in Eq. (^). 
The correct size of /i-term is obtained if A ~ 10^^ GeV. The form of the Kahler potential 
is implicitly suggesting that the Higgs fields have some interactions with the supersymmetry 
breaking sector mediated by particles with masses of 0(A). We also assumed that there is an 
approximate PQ symmetry discussed in subsection |2.2| with PQ{S) = 2. With the PQ symmetry, 
we cannot write down any other terms. Since S carries a charge, the Kahler potential for S 
is restricted to be the form in Eq. (^). The term in Eq. (§) represents the explicit but small 
breaking of the PQ symmetry. 

We here summarize the set-up. We consider the effective Lagrangian written in terms of the 
Goldstino multiplet S and the MSSM matter/gauge fields: 

A2 



A / A2 



+ (l-^^pC2(i?)(log|5|)2)ci>t$ , 

(27) 



W = WYukawa(^) + S + Wq , 
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The chiral superfield $ represents the matter and the Higgs superfields in the MSSM, and 
W^Yukawa IS the Yukawa interaction terms among them. We defined 0(1) valued coefficients 
cs, c^, and ch- We normahze the A parameter so that C5 = 1 in the following discussion. 
The parameters ch and A take real values whereas is a complex parameter. We consider 
the supergravity Lagrangian defined by the above Kahler potential K, superpotential W, and 
gauge kinetic function /. This is a closed well-defined system. The linear term of S in the 
superpotential represents the source term for the i^'-component of S. The last term in the 
superpotential, wq, is a constant, \wo\ ~ m^Mpi/\/3, which is needed to cancel the cosmological 
constant. The scalar potential has a minimum at (S) ~ A^/Mpi which avoids the singularity 
at S = 0. The set-up includes the dynamics of supersymmetry breaking and mediation. By 
expanding fields from their vacuum expectation values, we can obtain all the mass spectrum 
and interaction terms. 

When we write down the Lagrangian of the standard model we usually include the Higgs 
potential, V = (A///4)(|ffp — f^)^, and the gauge interaction terms of the Higgs boson instead 
of just giving bare mass terms to the W and Z bosons. Analogous to that, the system above 
contains dynamics of the supersymmetry breaking and a mechanism of its mediation instead of 
simply writing down soft supersymmetry breaking terms.f] In this sense, this way of construction 
is essential for the model to be called the MSSM in a true meaning. 

The effective Lagrangian is defined at the scale where the messenger fields are integrated out. 
The messenger scale, k{S), is not necessary to be 0{{S)). The k parameter originally comes 
from superpotential terms like, W 3 kSff. If the S field is a composite operator above the 
scale A as is often the case in dynamical supersymmetry breaking scenarios, the k parameter is 
suppressed by a factor of (A/Mpi)'^^'^)"^, where d{S) is the dimension of the operator S above 
the scale A. Therefore, the size of k depends on the actual mechanism of the supersymmetry 
breaking. 



We can see very nontrivial consistencies in this simple set-up. First, the /x-term is generated 
by the Kahler term, HuHd/A: 



''Our construction should not be confused with the spurion method of writing down the soft terms. The field 
5 is a propagating field and obeys the equation of motion. 




(28) 



With the shift of {S) in Eq. (|2^), the gaugino masses are 




(29) 
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Figure 2: Phenomenologically required values of the Higgsino mass ]x (with an 0(1) ambiguity, 
see text), the Bino mass and the gravitino energy density 03/2/1^- These three quantities have 
different dependencies on parameters 777,3/2 and A. The three bands meet around 7713/2 ~ 1 GeV 
and A ~ Mqut- The quantity VL^i2h? is defined in Eq. (30). It represents the energy density 
of the non-thermally produced gravitinos through the decays oi S \i S ^ hh is the dominant 
decay channel. 



Here and hereafter, we take a minimal model with A^n^ess = 1- The qualitative discussion does 
not change for different values of A'mess- Similar sizes of scalar masses are obtained from the 
Kahler terms. Finally, the moduli problem now turns into a mechanism for the production of 
dark matter. The energy density of the coherent oscillation of S dominates over the universe, 
and the reheating process by decays of the S'-condensation later produces gravitinos through a 



rare decay process S — > 7/^3/2 V'3/2- The amount can be expressed in terms of 7773/2 and A [17|: 

,3/2 / A ^ 



^3/2/1^ = 0.1 X 



777 



3/2 \ 



500 MeV 



1 X 1016 GeV 



(30) 



Here we have assumed that the decay of S into two Higgs bosons, S hh, is the dominant decay 
channel. The phenomenological requirements that /j, ~ 777^^/2 ~ O(IOO) GeV, and $13/2/1^ — 0.1 



can all be satisfied when 777 



3/2 



1 GeV and A ~ 10^^ GeV. 



We can see the non-trivial success of this framework in Fig. ||, where we see how 0(1) GeV 
gravitino mass is selected. The bands of 100 GeV < p, < 500 GeV, 100 GeV < 777^ < 500 GeV, 
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and 0.08 < ^^3/2/1^ < 0.12 are shown, where we defined /Z = mg/sMpi/A and ^3/2/1^ by Eq. 

The Bino mass is the mass of the U(l)y gaugino. Surprisingly, these three bands meet at 

m3/2 ~ 1 GeV and A ~ Mqut ~ 10^*^ GeV. 

The fact that A coincides with the unification scale, Mqut, is also quite interesting. In 
grand unified theories (GUTs), such as in SU(5) or SO(IO) models, we need to introduce colored 
Higgs fields in order for models to be consistent with gauge invariance. The colored Higgs 
fields, however, need to get masses through the spontaneous breaking of SU(5) or SO(IO). This 
suggests that the Higgs multiplets have some interactions with the GUT-breaking sector whose 
typical mass scale is, of course, Mqut- Therefore, it is quite natural to have Mqut suppressed 
interactions in the low energy effective theory. The same "cut-off" scale A for S suggests that the 
dynamics of GUT breaking is responsible for the supersymmetry breaking as well. The picture 
of unification of the Higgs sector, the supersymmetry breaking sector and the GUT breaking 
sector naturally comes out. Although it sounds like a very ambitious attempt to build a realistic 
model to realize this situation, it is quite possible and even very simple to build such a dream 
model by using a recent theoretical development of supersymmetric field theories [^6|. For an 
explicit example of such a GUT model, see Ref. [^] . 

This quite simple framework summarized in Eq. (pTj), gauge mediation with direct couplings 
between supersymmetry breaking sector S and the Higgs fields at the GUT scale, solves all the 
problems we mentioned before. We discuss these one by one here. 



Supersymmetric flavor problem 

The gravity mediated contributions to the sfermion masses squared are of O {171^2) ■ Therefore, 
the flavor mixing in sfermions are at most of 0(10^'*) level for 171^/2 ~ 1 GeV. For example, the 
constraints from the fj, ^ ej decay and /x — > e conversion process in nuclei put bounds on the 



mixing to be |12] 



(4)iV~f^)(^=^);sio-«, (31) 

V"^SUSY / V rnsvSY ) 

where msuSY is a typical sfermion/gaugino mass scale and is the muon mass. By using 
the fact that the value of tan/3 (= {Hu) / {Hd)) is predicted to be 0(30 — 40) (see discussion 
in the next section), this bound is marginally satisfied with sfermion masses of O(IOO) GeV. If 
gravitational dynamics maximally violates flavor conservation, future or on-going experiments 



have good chances to see the effects |35 



"The band of $73/2^'^ does not represent the dark matter density once we deviate far from the region of 
n^3/2 ~ 1 GeV. For rrig^p, or rrig^Jj,, the successful electroweak symmetry breaking cannot be achieved, and 
we cannot perform a sensible calculation of the 5 — > fefe decay width. 
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The flavor mixings from the high-scale dynamics such as physics at the GUT scale and 
the effect of right-handed neutrinos [^] are small as is always the case in gauge mediation. 



Super symmetric CP problem 

There are two physical phases in the MSSM: 

arg(mi/2/i(S/i)*) , arg(mi/2^*) . (32) 



From the Kahler term in Eq. (|26D, A- and i?-terms of 0(7713/2) are generated, but these will be 
overwhelmed by one-loop renormalization group (RG) contributions below the messenger scale. 
Since the RG contributions are proportional to the gaugino masses, the physical phases above 
are approximately vanishing. 



In fact, the phases of the original 0(7/73/2) contributions are also aligned with those of gaugino 
masses. The phases of three complex parameters in the Lagrangian, i.e., m?, c^, and wq, can 
all be taken to be the same by a field redefinition via and U(1)pq transformations. 

Even if there are 0(1) phases in the 0(7)7,3/2) contributions, which is possible if the PQ 
symmetry is maximally violated by operators suppressed by the Planck scale, the above physical 
phases are of 0(1%) which again marginally satisfies the experimental constraints. The upper 



bound on the electric dipole moment of the electron, for example, gives a constraint [12|: 

' rrie tan /? 



/ "73/2 \ 



< 10-^ (33) 



V"^SUSY/ V "T-SUSY 

where 777e is the electron mass. The bound corresponds to 77isusY ^ 300 GeV for 7773/2 ~ 1 GeV. 
/i-problem 

There are three kinds of //-problem in the MSSM, i.e., " Why fi < Mpi?", " Why fi^ ~ m^J", 
and " Why fi ~ 7771/2?" The second and third ones are related because there is a one-loop 
correction to the ni'jj^ parameter proportional to "^1^2' 

The first one was answered by the approximate PQ symmetry. The /7-term is forbidden by 
symmetry, but induced by a small explicit breaking term, W B m'^S. 

We can naturally obtain the relation, fi'^ ~ th,^^, once we assume the form of the Kahler 
potential to be the one in Eq. (|26|). This is a generalization of the Giudice-Masiero mechanism 



in gravity mediation [11|. The relation is independent of the "cut-off" scale A. We discuss a 



possible origin of the Kahler terms later. 
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The final relation, /i ~ is realized when A ~ Mqut as we can see in Fig. ^. Prom 

Eq. ( p8[ ) and (|29|), the relation between p. and the Bino mass, m^, is 

= 0-6 X fTTTTT^TT-TTl • (34) 



^lxl0i6GeV^ 

Although it is an 'accident' to have similar values of fi and the gaugino masses, the value we 
need, A ~ MguTj is motivated by two other independent physics, i.e., grand unification and 
dark matter of the universe. 

Cosmological moduli/gravitino problem 

As we have already discussed, the energy density carried by the coherent oscillation of S would 
not cause a problem. The decay of S reheats the temperature of the universe to of O(IOO) MeV 
for ~ 1 GeV and A ~ 10^^ GeV. This is high enough for the standard BBN. The non- 

thermal gravitino production from this decay gives the largest contribution to the matter energy 
density of the universe. The amount in Eq. (^) is, amazingly, consistent with the observation. 

The baryon asymmetry existed before S decays is diluted by the entropy production. If we 
assume the initial amplitude of S to be of 0(A), the dilution factor is estimated to be of order 
10~^(Tfj,/10^ Qgy-j-i ^[^]^ j'j^ ^j-^g reheating temperature after inflation. Therefore, a larger 
amount of baryon asymmetry is needed to be generated if baryogenesis happened above the 
temperature of O(IOO) MeV. 

If the stau is the NLSP as in the case we will study later, staus are also non-thermally 
produced through the S decays if it is kinematically allowed. The pair annihilation process 
reduces the amount but the abundance ends up with of 0(50) times larger than the result of 
the standard calculation of the thermal relic abundance. There are constraints on the decay of 
the staus into gravitinos from the BBN. Recent calculations including the catalyzing effects give 
an upper bound on the life-time of stau to be 0(1000) seconds |42, 43, 44, EH]. Although the 



lifetime is extremely sensitive to the stau mass (oc the typical lifetime with ~ 1 GeV is 
on the border of this constraint. This coincidence may be interesting for the Lithium abundance 
of the universe MS, Eql . 



Unwanted axion? 

It is common in gauge mediation scenarios that there is an approximate U(l) symmetry which 
is spontaneously broken. Therefore there is a (possibly unwanted) Goldstone boson associated 
with it |l^]. In the scenario we are discussing the vacuum expectation value of {S) breaks the 
approximate PQ symmetry spontaneously. The axion associated with the symmetry breaking 
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is actually the scalar component of S itself. The S scalar has a mass of the order of 100 GeV 
(see Eq. (^)) because of the linear term in the superpotential. Interactions between PQ currents 
and the axion S are suppressed by the scale of the symmetry breaking {S) ~ 10^^ GeV. There 
is no experimental or astrophysical constraint on such a particle. As we discussed above, the S 
scalar even plays an essential role in cosmology. 

Dimension-four and five proton decay problem 

The dimension-four operators which violate the baryon number conservation are forbidden by 
an unbroken Z2 subgroup of the PQ symmetry. This is identical to the i?-parity. 

Dimension five operators, such as QQQL, are allowed to appear at low energy because 
the PQ symmetry is spontaneously broken. In particular, if there are following terms in the 
superpotential: 

SQQQL , SUUDE , (35) 

the dangerous terms like QQQL and UUDE appear by substituting the vacuum expectation 
value of 5" ~ A^/Mpi ~ 10^^ GeV. In GUT models, these effective operators can be generated 
by diagrams with colored-Higgs exchange. In this case, the coefficients of the above operators 
will typically be of 0{fufd/ ^gut) where /„ and fd are the Yukawa coupling constants of 
up- and down-type quarks. By substituting (5), this becomes effectively QQQL or UUDE 
operators suppressed by fufd/Mpi. The prediction to the proton life-time is on the border of 
the experimental constraints with such coefficients . 



UV completion 



The discussion so far is based on the low energy effective theory defined in Eq. (27). This 
effective theory is valid up to the messenger scale k{S). Although it is not necessary for the 
discussion of low energy physics to specify UV models, an existence proof of an explicit UV 
completion supports our ansatz in Eq. (^). 

It is straightforward to UV complete the theory above the messenger scale by simply assuming 
a presence of messenger particles / and / which carry the standard model quantum numbers, 
and an interaction term kSff. The full model is K B f^f + p f and W 3 kSff instead of 
terms involving log S in Eq. 



The model with messenger fields now has a supersymmetric and hence stable vacuum at 



S = and f = f = \l —m?/k. However, as it has been shown in Ref. |22], there is a meta-stable 
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Figure 3: Feynman diagrams to generate higher dimensional operators in a UV model. 

minimum at (S) ~ A^/Mpi where supersymmetry is broken and messenger fields are massive. 
The effective theory in Eq. (^) correctly describes physics around the meta-stable vacuum. 

Above the mass scale A, we need a further UV completion. The simplest model is the 
O'Raifeartaigh model ||Tot| : 

K = S^S + X^X + y V , (36) 

and 

Ws = m^S + "^SX^ + MxyXY , (37) 

where k and Mxy{^ ^) a coupling constant and a mass for X and Y, respectively. There 
is an approximate PQ symmetry with charges PQ{X) = —1 and PQ{Y) = 1. By integrating 
out massive fields X and Y, we obtain the Kahler term — (S'^5)^/A^ with 

' "^'^ ' (38) 



A2 12(47r)2M2.y ' 

at one-loop level (see Fig. ^). The Higgs fields can directly couple to this system so that we obtain 
effective operators in Eq. (^). The terms are generated by introducing following interaction 
terms in the superpotential: 

Wuiggs = hHuqX + hHaqX + M^qq , (39) 

where h and h are coupling constants. Again, the PQ symmetry is preserved for PQ{q) = 
PQ{q) = 0. The supersymmetry breaking still happens in this extended model. After integrating 
out q and q, we obtain the Cf^S^ HuHd/ A term with 

A (47r)2 M, ■ n ; ' ^^^^ 



where 
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The term is also generated with 



A2 (47r)2 M2 ^ V ^1 
where 

-3 + 4x-x2 -2Ioga; 
^^(-) = ^(13^)3 • (43) 

These are obtained by calculating Feynman diagrams in Fig. |^. No other unwanted terms are 
generated because of the approximate PQ symmetry in the model. 

We can obtain the relation ~ c// ~ 1 for Mxy ~ Mq if the values of k, h and h are 
relatively large. In particular, we find 

where x = M'^y/^q- The ^-term squared is suppressed by a one-loop factor compared to the 
soft mass term 'm'jj^ for Mxy ~ Mg. The function f{x)^/g{x) never exceeds 0(1) values even 
for general relations between Mxy and Alg. To avoid a too large hierarchy, the loop expansion 
parameter (47r)2 should not be too small, i.e., the model should be (semi) strongly coupled.^ 
This fact suggests that this O'Raifeartaigh model itself is an effective theory of some dynamical 
supersymmetry breaking models. 

Indeed, there is an incredibly simple dynamical model which provides the above O'Raifeartaigh 
model as an effective description. The model is also embeddable into an SU(5) unified model in 
a straightforward way. The same dynamics spontaneously breaks SU(5) gauge symmetry and 
supersymmetry [37|. 



The model is based on a strongly coupled gauge theory where S and the Higgs fields appears 
at low energy as massless hadrons. The constituent 'quarks' of these hadrons are Q, Q and T, 
all of which transform as a vector representation of a strong S0(9) gauge group. Also Q and Q 
carry standard model quantum numbers (5 and 5 under SU(5)) and T is singlet under SU(5) 
but carries PQ{T) = 1. (See Fig. |^ for the structure of the model.) The Higgs fields and S are 
identified with meson fields 

H ~ {QT) , H ~ (QT) , S ~ (TT) , (45) 

where the Higgs fields in the 5 and 5 representation, H and H, contains and as SU(2) 
doublet components, respectively. 



"Note that this is not the same situation as the discussion around Eq. (20). The one-loop factor enhancement 

1/2 

9- 



there, Bfi/fi — mi/2/(ff^/(47r)'^), is always large due to the perturbativity of the standard model gauge coupling 
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S0(9) becomes strong 
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Figure 4: Structure of an example of the UV model 

This is an SO (9) gauge theory with eleven flavors, and the SU(5) gauge group is identified 
with a subgroup of the SU(ll) flavor symmetry. We can write down superpotential terms: 

T^GUT = MtT^ + MqQQ - (QQf + ■■■ , (46) 

where //^ (~ 1 — 10 GeV) corresponds to the small explicit breaking of the PQ symmetry. This 
term is going to be the m'^S term in Eq. (|27| ) at low energy. Once we ignore the superpotential 
(in the limit of fix^Mq and Mx oo), the S0(9) 11 flavor theory is on the edge of 
the conformal window [^8|. Therefore, at some scale A* the gauge coupling constant flows 
into the infrared fixed point. Although it becomes a strongly coupled conformal field theory 
(CFT) near the fixed point, there is a dual weakly coupled CFT description with which we can 
perform perturbative calculations. The dual gauge group is S0(6) and the superpotential above 
is replaced with 

W^GUT = ^tA*5 + MqA,M-4^M2 + ... 

Mx 

+ y + KHqt + KHqt + KMqq + ■■■ (47) 

The field M is a composite meson M ~ {QQ) which transforms as 1 + 24 representation under 
SU(5). The fields t, q, and q are dual quarks which are charged under S0(6), and and /i* are 
the coupling constants at the fixed point (k,,, = /i* ~ (47r)/A^ with N = 6). At a vacuum where 
the gauge group is broken down to the standard model gauge group, (M) = diag.(0, 0, 0, w, "L") 
and {qc) ^ {qc '■ colored components of q), this model becomes exactly the same as the 
O'Raifeartaigh model in Eqs. (|37| ) and (^) with the identification oft ^ X, He Y {Hq ■ the 
colored Higgs field), /U^A* —>■ m?, h^{qc) Mxy and h^{M) Mq. (See Fig. ^ for particles 
to describe the effective theory in each energy interval.) 

Once we take into account non-perturbative effects, there appears a supersymmetric mini- 



mum far away from the origin of S. However, it has been shown in Ref. [26| that the vacuum 
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Figure 5: Particles to describe the theory in each energy interval. 

near S" = is meta-stable. We can see in Eq. ( p8| ) that the S mass squared, m| = +4F|/A^, is 
indeed positive. 

The loop expansion parameter, |/i*pA^/(47r)^ where = 6, is 1/Np at the fixed point in this 
model (Np = 11). Therefore, we obtain 

.2 1 

[Np = 11) . (48) 



m 



Hu 



1 

Np 



Although it looked problematic to have a hierarchy in Eq. (^) in perturbative models, similar 
sizes of /i and mn^ can be obtained in this semi strongly coupled theory: fi/mH^ ~ 1/3 

Doublet-Triplet splitting problem 



The model above completely solves the doublet-triplet splitting problem in GUT models. By the 
vacuum expectation value of the colored component of q and q, the gauge group SO (6) x SU(5) 
is broken down to the standard model gauge group. The Hqt and Hqt couplings in Eq. ( |47| ) 
then give mass terms only for the colored Higgs fields ||3^. This dual picture is similar to an 
SO(IO) model proposed in Ref. | |5]| . 

As discussed before, the dimension five operators for proton decays are sufficiently suppressed 
thanks to the approximate PQ symmetry. 



^^The relation is not a precise prediction of the model. Depending on the ratio of the mass parameters 
Mq(= h,{M)) and Mxy{= h^^qc)), which are independent parameters in the superpotential, there can be 0(1) 
deviation from the relation (see Eq. (^)). If Mq < Mxy , we can reliably use Eq. ^A) (by multiplying a factor of 
A'^) with h^N / [A-kY ~ 1/Nf as the leading order result of the 1/A'' expansion. However, once Mq/MxY becomes 
too large, such as a factor of three or so, we lose a perturbative control of the calculation. In this case, we can 
first integrate out qo and qo (the doublet part of q and q), and then by taking the Seiberg duality |49| of this 
dual picture again the theory becomes (semi) weakly coupled. (It is an S0(5) 7 flavor model. See Ref. [^.) 
Although there are 0(1) ambiguities in model parameters through the matching between two theories, the naive 
dimensional analysis fcol] gives the same result as the relation in Eq. (Uq) even in that case. 
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Supersymmetric fine-tuning problem 



The experimental lower limit on the Higgs boson mass from LEP-II experiment, > 114 GeV [52|, 
has put a threat on supersymmetric models. In order to satisfy the experimental bound, we 
need either a heavy scalar top quark (stop) or a large At-term (the stop-stop-Higgs coupling) 
since a significant one-loop contribution to ruh is necessary On the other hand, once we 

have large or At, it induces a large one- loop contribution to the soft mass term rn'jj^. This 
immediately means that there is a fine-tuning in the electroweak symmetry breaking as we can 
see in the condition: 



M2 

llLZ. ~ ..2 _2 f^^ 



2 - fi'-mi,^{A)-6mij^ , (49) 

where Srnjj^ is the contribution from the radiative correction, and Mz is the Z boson mass 
{Mz = 91.2 GeV). If |5m^^| ^ we need cancellation between Srnjj^ and either or 

m'jj^{A) to reproduce a correct value of the Z boson mass. A cancellation of at least 0(1 — 5%) 
is necessary to satisfy the bound on the Higgs boson mass in generic gravity or gauge mediation 
models (see for review [^). 



Although the framework in Eq. (27) does not avoid the problem, there is an interesting 
consistency. As we have observed in an example of the UV completion before, the ratio 
of / m'jj^{A) is predicted to be small (at least a factor of a few) if the theory has (semi) 
perturbative description above the scale A. In general, without specifying UV models there 
is a good reason to believe that the description should be (semi) perturbative. First, we are 
implicitly assuming that quarks and leptons remain to be elementary particles (weakly coupled) 
all the way up to the Planck scale, otherwise we reintroduce the flavor problem. If quarks and 
leptons are strongly coupled above the scale A, it is expected to have interaction terms such 
as S^^^/ A"^ which induce mixing terms of 0(1) in sfermion mass matrices. On the other 
hand, we must write down the Yukawa coupling constant for the top quark which is of 0(1). 
This indicates that the Higgs fields should not be replaced by a composite operator with a large 
dimension. If the dimension of the Higgs operator above the scale A was d{H) > 1 and the top 
quark is an elementary particle as discussed above, the Yukawa coupling would be suppressed 
by a factor of (A/Mpi)'^^^^~^ . The 0(1) Yukawa coupling constant suggests that dimension 
of the Higgs operator is d{H) ~ 1, i.e., the Higgs fields are not (very) strongly coupled. The 
loop expansion should then make sense in the calculation of fi and m|^^(A). The smallness of 
fi^ / m'jj^{A) is, therefore, a generic feature of the model. 

Another amusing point to notice is that the function g{x) in Eq. ( |4^ ) is positive valued for 
X > 0. This means m'jj (A) > if there is a (semi) perturbative description. 
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Now with positive m'jj^{A) and fj,^ <^ m^^(A), the condition of the electroweak symmetry 
breaking in Eq. ( |49D imphes that S^h^ must be negative and large. Indeed, the contributions 
from the stop-loop diagrams are negative and are proportional to m?. Therefore, the "little 
hierarchy", i.e., a heavy stop is predicted in this model. The tight experimental bound on the 
Higgs boson mass is not a big surprise. 



Strong CP problem 

Although the approximate PQ symmetry introduced in the framework does not provide us with 



a solution to the strong CP problem in a usual way by the axion mechanism [55| (because it is 
explicitly broken), there is an interesting connection. 

The PQ symmetry is anomalous with respect to the SU(3) strong interaction of the standard 
model. If we demand the PQ symmetry to be non-anomalous, there are two options to take. 
The first one is to simply assume that the u-quark is massless. By combining U(1)pq with the 
chiral symmetry, we can make the PQ symmetry non-anomalous. 

Another option is to introduce an axion chiral superfield A which has a coupling to the gauge 
fields: 

/ 9 4-W"Wa . (50) 
JA 

The kinetic term for A is 

K3{A + A^f. (51) 

With a PQ transformation A, A ^ A + i6, we can cancel the gauge anomaly. 

Each of two options, massless u-quark and the axion, solves the strong CP problem.^ 



3 Low energy predictions 



The set-up in Eq. (p7|) provides a characteristic spectrum of the supersymmetric particles. It is 

different from conventional gauge or gravity mediation models. Since the Higgs sector directly 

couples to the supersymmetry breaking sector at the GUT scale, the soft mass terms for the 

Higgs fields are generated at the GUT scale. The gaugino masses and sfermion masses are, 

on the other hand, generated at the messenger scale. This hybrid feature provides interesting 

predictions on the low energy spectrum. 

*The scalar component of the axion chiral superfield A, the saxion, obtains a mass of the order of 1 GeV by 
gravity mediation effects |^^. With a sufficiently small decay constant /a, the saxion does not cause a moduli 
problem as it decays into a pair of gluon much earlier than the decay of S. 
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We discuss a parametrization of the model defined in Eq. (^T\), with which we can calculate 
the low energy spectrum and interaction terms. As we will see below, we can parametrize the 
model by three quantities. These three define a theoretically well- motivated hypersurface in the 
large dimensional MSSM parameter space. 

3.1 Parametrization 

We first count the number of the parameters in the model. The soft supersymmetry breaking 
terms for the Higgs sector: 

mjj, n , (52) 

are generated at the scale A. We take the scale A to be the unification scale Mqut- We assumed 
the same soft mass terms for Hu and Hd (m^^(MGUT) = "^^^(-^gut) = itt-h) as motivated by 
the UV completion discussed before. Gaugino masses, ^-terms, i?-term, and sfermion masses 
are vanishing at the GUT scale. 

Below the GUT scale, RG evolutions of the soft terms induce sfermion masses through the 
Yukawa interactions. The gaugino masses. A- and i?-terms remain vanishing. At the messenger 
scale, 

Mmess , (53) 

the messenger fields decouple. The threshold corrections (i.e., gauge mediation effects) con- 
tribute to the gaugino masses, sfermion masses and also the Higgs masses squared. Those are 
calculable with a single parameter. 



as we can read off from Eq. (27). This parameter controls the overall scale of the supersymmetry 
breaking parameters. For example, the gluino mass is M3 = c/|M gl, |. The A- and B-terms 
are still vanishing (up to higher order loop corrections IS^) at the messenger scale, but the RG 



evolution below the messenger scale generates those through one-loop diagrams. 

All the soft supersymmetry breaking parameters at the electroweak scale can be expressed in 
terms of these four parameters, m^, fx, Mmess, and M, by the procedure described above. One 
combination of the parameters should be fixed by the condition for the electroweak symmetry 
breaking, i.e., Mz = 91.2 GeV. We take the rnjj parameter as an output of the calculation. The 
model parameters are now defined by (^, Mmess, -^)- Here we take the running /x parameter at 
the scale MgusY = ("if w-f )^^^ as an input parameter. 
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Figure 6: The RG evolution of the super symmetry breaking parameters. RG equations at one- 
loop level are used. A parameter set (/x, Mmess,M) = (300, 10"'^'^, 900) [GeV] is chosen. The 
left panel shows the evolution of the soft masses for iji (dotted), fji (dot-dashed), and Hu 
(solid). The rhx parameter is defined by rhx = sgn(m^)|m^|^/^ for each chiral superfield 
X. The evolution of the /^-parameter (dashed) is also shown. Negative contributions for fh^ 
and fhf above the messenger scale comes from the one-loop contribution through the Yukawa 
interactions. Threshold effects (gauge mediation) at the messenger scale contribute to sfermion 
and the Higgs mass parameters. The ni'jj^ parameter is driven to a negative value by the 
stop-loop diagrams. In the right panel, gaugino masses, A-, and S-parameter are shown. The 
gaugino masses are generated at the messenger scale, and induces A- and i?-terms by the one- 
loop running. For the phase convention of A- and i?-terms, we have used the one defined in 



Ref. |8|. 



We show in Fig. ^ an example of the RG evolution of soft supersymmetry breaking parameters 
for (/i,Mmess,M) = (300 GeV, 10^° GeV, 900 GeV). The horizontal axis hr is the RG scale. 
We have used the top quark mass, rrit = 170.9 GeV |59| |. The constraint from the electroweak 
symmetry breaking fixes the m'jj parameter to be (817 GeV)^. The choice of parameters is 

^1 



motivated by the discussion in the last section. The positive value of m|r and a relatively small 



value of //(Mgut) compared to \Jm?^ are realized with this set of parameters. The lightest 
Higgs boson mass is calculated to be 115 GeV. We will use this set of parameters in a collider 
study in Section 

In the left panel of Fig. ^, scalar masses and the /^-parameter are plotted. We have defined 
mass parameters rhx = sgn(m^)|m^|^''^ for each scalar mass parameter m?x- Several interesting 
things are happening here. With non-zero positive values of m?^ the Yukawa interactions induces 
negative masses squared for sfermions in the third generation. The positive values are motivated 
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by the positivity of the function g{x) in Eq. (^3|). The negative contributions to the sfermion 
masses are compensated by the positive contributions from the gauge mediation effects at the 
messenger scale. This behavior of the RG evolution gives smaller values of the stau mass, ^f^, 
at the electroweak scale compared to those in the conventional gauge mediation scenarios. The 
impact on the collider physics of this effect will be discussed in Section 0. 



As is clear from Eq. ( p9D the value of the //-parameter is approximately obtained by /x ~ 
—'f^Hu Ufl^TeV- Therefore, we can easily understand from this figure that the value of /j, is smaller 
for larger initial values of rnjj. On the other hand, the stau mass is also smaller for large m'jj 
because it receives more negative contributions. This correlation between the /^-parameter (the 
Higgsino masses) and the stau masses is an interesting prediction of the model. We will discuss 



this relation more in subsection |3.3 . 

The evolution of gaugino masses, A- and S-terms are shown in the right panel of Fig. ^. 
We have used the sign convention of those parameters defined in Ref. [^]. Those i?-charged 
parameters are generated at the messenger scale by gauge mediation effects. We see a peculiar 
behavior of the B-term. It starts from zero at the messenger scale, goes to negative once and flips 
its sign to positive later. The negative contribution comes from a loop diagram with the SU(2) 
gaugino, and a stop-loop diagram with the At term gives a positive contribution. The behavior 
can be understood by the fact that the At term starts at zero but its absolute value becomes 
larger at low energy. The positivity of the i?-term remains to be true unless the messenger scale 
is extremely low such as Mmess ^ 10^ GeV. In this case, the two physical signs are predicted to 
be 

sgn(Mi//(B//)*) = +1 , sgn(Mi^*) = -1 {i = 1, 2, 3) . (55) 

Unlike other models used in literatures for collider studies, there is no choice of these signs. The 
former predicts a positive sign for the supersymmetric contributions to the anomalous magnetic 
moment of muon. Also, the latter determines the sign of the chargino-loop contribution to the 
amplitude of the 6 ^ 57 decay to be opposite to the one from loop diagrams with the charged 
Higgs boson. Interestingly, both of these signs are preferred by the experimental constraints on 
these quantities. 



3.2 Electroweak symmetry breaking 

It is highly non-trivial whether we can have successful electroweak symmetry breaking with the 
limited number of parameters. We demonstrate here that the correct Z-boson mass can be 
obtained without spoiling the perturbativity of the Yukawa interactions for the top and bottom 
quarks up to the GUT scale. 
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H [GeV] H [GeV] 

Figure 7: The corresponding value of rhn (left) and tan/3 (right) to the input parameter 
(/x(Msusy))- We set the overall scale M = 900 GeV. For different values of M we can obtain 
approximate relations by rescaling the axes. Curves for messenger scales Mmess = 10^, 10*^, 
lO^'', and 10^^ GeV are shown. The curves are terminated by the mass bound of stable staus 
nif-^ > 98 GeV [^]. Small values of ^'^/mjj are predicted if the UV theory is weakly coupled. A 
rough classification of 'weakly coupled', 'semi perturbative', and 'strongly coupled' is indicated. 

The left panel in Fig. |^ shows the value of friH (defined at the GUT scale) required by 
the correct electroweak symmetry breaking with respect to the running ^u-parameter at MsusY- 
Curves for different messenger scales Mmess are shown. The M parameter is fixed to be 900 GeV. 
For other values of M, say xM with an arbitrary positive number x, we can obtain curves, as a 
good approximation, by rescaling the both axes by the factor of x. As is clear from the behavior 
of the RG evolution shown in Fig. ^, larger values of fhu are necessary for having smaller values 
of fi. Each lines are terminated at some small value of /u, where the stau becomes too light 
(mfj < 98 GeV [|6^) due to the negative contribution from the running between the GUT scale 
and the messenger scale. This negative contribution is significant only when tan /? is large. It is 
indeed the case as we will see later. 

As we discussed in the previous section, the UV completion of the theory above the GUT 
scale suggests that the Higgs fields do not get strongly coupled. The discussion is based on the 
requirement of the 0(1) Yukawa coupling constant for the top quark. This indicates a (small) 
hierarchy between the m//-parameter and the /u-parameter at the GUT scale since the ratio 
fi^/mjj turns out to be the loop-expansion parameter. We indicate in the plot the region with 
< fi^/m% < 1/100 to be 'weakly coupled', 1/100 < fi^/mjj < 1/4 to be ' semi perturbative', 
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and fi'^/m'jj > 1/4 or < to be 'strongly coupled' for illustration. (There is no concrete 
meaning in precise locations of border lines. )[| We only find solutions with small fi compared to 
the overall scale M if we impose (semi) perturbativity. The light Higgsino is, therefore, one of 
the predictions of the model. Also, a high messenger scale, such as Mmess ^ 10^ GeV is required 
to be in that region. 

There is a possibility that the Higgs fields are fully strongly coupled and the top Yukawa 
coupling is generated by some strong dynamics by making the top quark also involved in the 
strong sector. Since there is no stringent constraint for flavor violation in the top-quark physics, 
this is not a dangerous assumption although the embedding to GUT models may be more 
difficult. In this case, the stop mass squared at the GUT scale becomes an additional parameter 
of the model. We do not consider this modification of the framework in this paper. 

We show in the right panel the predicted values of tan/3. Again, we take M = 900 GeV. 
The rescaling of the horizontal axis gives a curve for different values of M with a good accuracy. 
For small values of /x preferred by UV completions, tan /3 ~ 30 — 40 is predicted. The values 
in this range do not cause a Landau pole of the Yukawa coupling constants below the GUT 
scale. Therefore, the perfectly consistent electroweak symmetry breaking is achieved without any 
extension of the model. Large values of tan/3 are due to the fact that the Bfj, term is generated 



only through radiative corrections below the messenger scale |57]. The Yukawa coupling constant 



for the tau lepton is large for large tan (3, which affects the running of m? . 
3.3 Light f and light Higgsino 

The characteristic RG evolution of the supersymmetry breaking parameters in this framework 
provides interesting predictions in low energy physics. In particular, the correlation of the 
^-parameter and the stau mass gives a large impact on collider physics. 

The plot in Fig. ^ shows the correlation for different values of Mmess- For small values of 
/i, the stau is lighter than the Bino and Higgsinos. Therefore, for values of /U in the (semi) 
perturbative region (see left panel of Fig. |^) there is a big chance that the stau becomes the 
NLSP (remember that the gravitino is the LSP). The collider physics of this region will be very 
different from the scenario with the neutralino NLSP as there is no missing Ej- associated with 
escaping neutralinos. (The lifetime of the stau is typically of 0(1000) seconds with which the 
stau is regarded as a completely stable particle for the time scale of collider experiments.) Also, 
the light Higgsino changes the pattern of cascade decays of supersymmetric particles. We will 



* Precisely speaking, the hierarchy is predicted for the ratio of ffiH and the jj, parameter at the GUT scale 
whereas what is plotted is the /ijarameter at Msusy- However, the RG running of the fj, parameter changes its 
value only by 0(10%) (see Fig. 
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Figure 8: Correlation of the fi parameter and the hghter stau mass. The overall scale M = 
900 GeV is set. Below the lines of the Bino mass and the Higgsino mass, the stau is the NLSP. 
We can clearly see the positive correlation. 

discuss in the next section the overall feature of this scenario at the LHC and demonstrate a 
way of confirming/excluding the framework. 



4 LHC signatures 



The theoretical success of the sweet spot supersymmetry motivates us to consider what will be 
the experimental signatures at the LHC experiments. We show in this section that there are 
several unique features. We present a way of confirming/excluding the model in the case where 
the lighter stau is the NLSP. 

4.1 Overview of supersymmetric events with r NLSP 

As we have seen in the last section, it is plausible that the lighter stau is the NLSP. A small value 
of the /x-parameter is a natural consequence of UV physics, and that makes f light through the 
RG evolution. If it is the NLSP, the lifetime of stau is of 0(1000) seconds with our assumption of 
the 0(1) GeV gravitinos. The LHC signals with such a long-lived stau will be quite different from 
ones with the usual assumption of the neutralino LSP. There have been many studies on collider 



signatures for the quasi-stable f-NLSP scenario, for example, in |61|-[72|. We demonstrate here 
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Table 1: Masses of superparticles and Higgs bosons in GeV for our benchmark point, /i = 
300 GeV, Mmess = 10^° GeV and M = 900 GeV. The gravitino mass is fixed to account for the 
observed dark matter density (see Eq. (|30|)). Here, the masses of the squarks and sleptons of the 
second generation are omitted, since they are equal to the ones of the first generation. We use 
the notation for the superparticles and Higgs bosons in the MSSM in Ref. |78]. 



reconstruction of model parameters with f NLSP at the LHC experiments. [| 
We select the following benchmark point for the collider study: 

M = 300 GeV , Mmess = 10^° GeV , M = 900 GeV . (56) 

This set represents the most theoretically motivated region of the parameter space as we have 
discussed before. As we can see in Fig. ^ the NLSP is the stau with this set of parameters. 
We have calculated the spectrum by solving RG equations at one-loop level. The running 
parameters at the scale MgusY = ("if "if = 1053 GeV have been used for the calculation 
of the spectrum. We have ignored the QCD finite corrections at the low energy threshold, 
which would amount to about 10%. In Table 0, we listed masses of superparticles and Higgs 
bosons. We used mt = 170.9 GeV. The stau mass is 116 GeV, and its lifetime is calculated 
to be 3000 seconds with the gravitino mass determined by the dark matter density, 771,3/2 = 
500 MeV. The running gaugino mass parameters at MgusYj Afj = gfM, are Mi = 195 GeV, 
M2 = 364 GeV, and M3 = 1013 GeV. The lightest neutralino Xi is, therefore, mostly the Bino, 
X2 and X3 mainly consist of the Higgsino components, and the Wino is the heaviest, Xi- The 
lighter and heavier charginos, Xi and xt^j are mainly the Higgsino and the Wino, respectively. 
The gluino and squark masses are about ITeV. The Higgs boson mass, 115 GeV, is calculated 



* Recent studies on the lifetime measurement of the long-lived charged NLSP in the collider experiments show 
that it is possible to determine the gravitino mass in some range of the parameter region although those 

proposals require an extra experimental set-up to collect the charged NLSP. 
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Figure 9: The dominant decay modes of g and q. The percentages show the branching ratios 
of each modes. The shaded modes are relevant for the analysis of the reconstruction of the 
neutralino masses. 




Figure 10: The dominant decay modes of x X ■ The shaded modes are relevant for the 
analysis of the reconstruction of the neutralino masses. 

by using a one-loop effective potential with taking into account leading two-loop corrections by 
appropriately choosing a renormalization scale for the running top quark mass which appears 
in the effective potential [^9|. Similar values, 114 — 115 GeV, are obtained by using publicly 
available codes ||8^, ^ . 

The total cross section of the superparticle production at the benchmark point is 1.4 pb for 
the center-of-mass energy of the LHC. The cross section is dominated by pair productions of 
gg, qg, and qq. The subsequent decays of these colored particles generate hard jets and other 
supersymmetric particles such as neutralinos and charginos. The decays of these non-colored 
superparticles, in the end, produce two quasi-stable fi's for each supersymmetric event. Most 
of the stau pairs escape a detector and leave two charged tracks. 

The decay cascades start with the the decays of g and q as shown in Fig. We have used 
ISAJET 7.69 to calculate the branching ratios. Since the gluino is lighter than squarks, it 
decays into a neutralino or a chargino through three-body decay modes. The dominant channel 
is the decay into a pair of third generation quarks and a Higgsino, Xi X231 through the 
Yukawa interaction of the top quark. The main decay mode of the squarks are q — > g + q, 
followed by the gluino decay. Therefore, for each supersymmetric event, many hard jets are 
produced. Especially, a significant number of 6-jets are produced by the gluino decays (and also 
by the subsequent decays of the top quarks). This is an interesting feature of the model, but 
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at the same time, the large number of jets makes it difficult to trace back the decay chains in 
the actual analysis because it is hard to specify which jet is the one originated from a particular 
decay process. 

At the end of the decay chain, three lighter neutralinos (Bino and Higgsinos) decay into fi 
and r (Fig. [l^) . The heaviest neutralino and the heavier chargino (Winos) decay into which 
in turn decays into a stau and a neutrino. We can measure the mass and the momentum of fi's 
in the final state by using information on the charged tracks in the muon system. With the full 
reconstruction of the four momentum of the staus, three of neutralino masses can in principle 
be measured as the invariant mass of fi and r. 

Sleptons except for fi do not appear in the decay cascades. That means we do not have 
clear lepton signals such as two opposite-sign same-flavor leptons from the X2 decay often used 



as a tool for precision measurements |83]. What we typically have are a lot of third generation 
quarks (6-jets) and leptons (r-jets) in the final states. It is again interesting but not exciting 
situation for measurements of the mass spectrum. 

4.2 Reconstruction of neutralino masses 

We demonstrate here a way of reconstructing the neutralino masses by using the decay modes 
X123 ~* '^i''"' ^'^d show that it is possible to determine the three model parameters. 

There have been studies on the mass reconstruction of the stau and the neutralinos in stau 
NLSP scenarios. Especially, in Ref. [^], a detailed study of the stau mass measurement is 
performed including detector effects. The method is measuring the velocity {(3) of the stau by 
the time-of-flight and the momentum {pr-^) from the track. We can then calculate the mass by 
the formula, 

mf, = ^- (57) 

It is concluded that the stau mass can be measured with the accuracy of at least about 100 MeV 
for TTif^ ~ 100 GeV. The standard model background is estimated in the paper. Since the 
staus with high velocities are indistinguishable from muons, a tight selection cut on the velocity, 
/37 < 2.2, is imposed in the analysis. We reproduced the analysis of the stau mass measurement 
and found that it is indeed possible to measure it with a very good precision at the benchmark 
point. We follow the selection cuts proposed in this paper in the analysis of the reconstruction 
of neutralino masses as well. 



The reconstruction of the neutralino masses has been discussed in Ref. |66| and also in a 



recent paper Q. Once we determine the stau mass, we know the four-momentum of the stau on 
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event-by-event basis. By combining with the four-momentum of r, we can extract the neutrahno 
masses. There are two things needs to be done in the analysis. Since every supersymmetric event 
contains two staus in the final state, there are two candidate staus to be combined with for each 
r. The invariant mass coincides with a neutralino mass only if we choose a correct combination. 
We need a strategy to select the correct one. The other thing is that we do not know the 
four-momentum of r. The r particle decays inside the detector and an invisible neutrino in the 
decay products always carries away some amount of energy. 



In the analysis in Ref. IjH^], it has been assumed that the selectron and the smuon are 
also quasi-stable. Therefore, by selecting events with only one stau there is no combinatorial 
background. A method of fully reconstructing the four-momentum of r has been discussed in the 
paper. If the direction of missing momentum is aligned with that of the r-jet, it can be identified 
with the neutrino. By adding missing momentum to the r-jet momentum, four momentum can 
be reconstructed. Even without trying to reconstruct the four-momentum of r, the authors 
found that there appear sharp edges at the neutralino masses in the distribution of the f — r 
invariant mass by using hadronic decay modes of r. 

We follow this endpoint analysis with hadronically decayed r's for the measurement of the 
neutralino masses. Since there are always two staus in the final state in contrast to the data set 
analyzed in Ref. [^], we need to consider a way of reducing the combinatorial background. 
Also, with updated experimental bound on the Higgs boson mass, the overall scale of the 
supersymmetric particles should be higher than the points studied in Ref. [^6| . This significantly 
reduces the statistics, which are essential for the endpoint analysis. Also, with the cut on the 
stau velocity discussed above, the number of signal events are again significantly reduced with 
the relatively light stau, 116 GeV, at the benchmark point. It is, therefore, non-trivial whether 
the method works in our model. 

It has been recently proposed to reconstruct the energy of neutrinos by decomposing the 
missing momentum into two directions of r's which we know from directions of leptons or r- 



jets [71|. Also, by using information of the electric charge of staus and leptons in the final state, 
we can select the correct combination for the events with f'*' and f~ . However, in our model, 
there are often other neutrinos in the events, e.g., from chargino decays and top-quark decays. 
The decomposition fails in the presence of such neutrinos. In the leptonic decays of r, neutrinos 
tend to carry away majority of r energy due to the kinematics and the V — A current structure 
of the weak interaction. Therefore, all of the uncertainties in the reconstruction of the neutrino 
momentum reflect to the mass measurement. For the background rejection, the authors have 
used a loose cut on the stau velocity, /37 < 6, and they checked that it is enough to reject the 
standard model background from mis-identifications of muons as staus. However, with a light 
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stau at our benchmark point, it is non-trivial whether such a loose cut can effectively suppress 
the background. The overlap between the signal and background regions in the (/37, mf) plane 
gets larger for a light stau. Once we impose a tighter cut, /?7 < 2.2, we could not obtain enough 
number of events for the analysis. Therefore, we do not pursue this direction. 

In the following, we discuss a method of the simulation and selection cuts for reducing the 



standard model background. The actual analysis will be presented in 4.2.2 
4.2.1 Event generation and selection cuts 

We have generated 42,900 events of the superparticle productions in proton-proton collision 
at the LHC energy by using a event generator HERWIG 6.50 Q with the MRST parton 
distribution function. The number of the events is equivalent to the integrated luminosity of 
30 fb"^ We have used the package TAUOLA 2.7 for r decays fsl. 



For a detector simulation, we have used a package AcerDET-1.0 |87], which is a fast simulator 



for high pt physics at the LHC. The AcerDET program identifies isolated leptons, isolated 
photons and isolated jets out of the final state of each generated event. The cluster selections 
and the smearing of the four-momentum of leptons, photons and jets are implemented. Leptons 
and photons are considered to be isolated if they are far from other clusters by AR > 0.4, and 
the transverse energies deposited in cells in a cone AR = 0.2 around the cluster are less than 
10 GeV. A cluster is recognized as a jet by a cone-based algorithm if it has pT > 15 GeV in a 
cone AR = 0.4. The package also implements a calibration of jet four-momenta using a flavor 
independent parametrization, optimized to give a proper scale for the di-jet decay of a light 
Higgs boson. Each jet is labeled either as a light jet, 6-jet, c-jet or r-jet, using information of 
the event generators. We have used default values of the parameters for clustering, selection, 
isolation, calibration, and labeling processes. For the r-jet identification, we further implement 
r-tagging efficiency of 50% per a r-labeled jet. 

A parametrization of the resolution of the stau velocity is shown in Ref. []6^ ]: 

^ = 2.8% X /3. (58) 

We smeared the velocity by using this resolution. The same resolution is assumed for the 
measurement of the muon velocity. 

For the smearing of the stau momentum, we have used the momentum resolutions shown in 



36 



Ref. 1 71 1]; one from the sagitta measurement error, 

0.0118% X (prJGeV), (59) 



Pfi 

one from a multiple scattering term, 




Pfi 

and one from the fluctuation of energy loss in the calorimeter. 



""^P^^^ -2%x Jl + -^, (60) 



criPfi] 



i X {prJGeY)-' (61) 



Pfi 

We have smeared the stau momentum according to these resolution width cr(pfi). 

If the measured velocity of the stau is high enough, such as (3j > 0.9 |7^, the stau will be 
identified with a muon and can be used as a trigger. However, for slow staus, we need to rely on 
other triggers. For the simulation of the triggering, we have chosen only events passing one of 



the following conditions [88|: one isolated electron with px > 20GeV, one isolated photon with 
Pt > 40GeV, two isolated electrons/photons with px > 15GeV, one muon with px > 20GeV, 
two muons with > 6 GeV, one isolated electron with px > 15 GeV and one isolated muon 
with Pt > Q GeV, one jet with px > 180 GeV, three jets with pT > 75 GeV, and four jets with 
Pt > 55 GeV. Here, isolated electrons/photons, isolated muons and jets must be in the central 
regions of pseudorapidity \r]\ < 2.5,2.4, and 3.2, respectively. Following |71|, we treated staus 
with > 0.9 as muons in the simulation of triggering. 

For the event selection, we require two stau candidates for each event. Since the stau mass 
can be precisely determined, a stau identification can be performed by testing if its measured 
mass by Eq. ( p7[ ) is consistent with the actual mass. For the consistency test, we took a window 
of the measured velocity, Pmeas- 

(3' - 0.05 < I3mc^, <P' + 0.05 , (62) 

where is a velocity calculated from the measured momentum, Pmeasj by assuming the stau 



mass, i.e., /?' = YPmeas/(Pmeas + "^ii) (sec |6^). To reduce the standard model background 
from mis-identifications of muons as staus, we required one of the stau candidate selected above 
to have < 2.2. The transverse momentum cut, pT > 20 GeV, is also imposed. The lower 
limit on the velocity /37 > 0.4 is imposed to ensure the stau to reach the muon chamber. As for 
the isolation of stau, we have used the same criterion with that of the muon. 



^According to the paper, the original study has been done by G. Polesello and A. Rimoldi, in ATLAS Internal 
Note ATL-MUON-99-06, but it is not publicly available. 
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The standard model background can be further suppressed by requiring the presence of the 
two charged stable particles as well as a large effective mass |^6|. The effective mass Mefj is 
defined by the scalar sum of the transverse momentum of the four leading jets pT,i (i = 1 — 4), 
and the missing transverse momentum E^^^^, 

Meff = PT,1 + PT,2 + PT,3 + PT,4 + ^t'"'- (63) 

The effective mass distribution of supersymmetric events has a peak around 1 TeV. We impose 
a cut Mcs > 800 GeV. With all of the selection cuts discussed above, the standard model 



background is reduced to a negligible level |68|. 



Events with only one r-tagged jet is selected for the reconstruction of the neutralino mass. 
We require pT > 40 GeV for the r-jet. If we allow two r-tagged jets, the number of signal 
events increases by about 10%, but it also increases the probabity of selecting a fake r-jet. The 
effects of the fake r-jets are taken into account by assuming the probability of the mis-tagging 
of non-r-labelled jets to be 1% per a jet |^^. We also discuss the significance of the fake events 
by comparing with a case with a 5% mis-tagging probability. For simplicity, we took the tau- 
tagging efficiency and the mis-tagging probability to be independent of jet px, rj, and the decay 
modes of r. 

4.2.2 Invariant mass analysis 

With the selection cuts described above, the total number of events are reduced to 2,000 in 
which we have 1,563 events with a true r-jet and a true stau pair. With these limited statistics, 
we need to develop an effective method to reduce the combinatorial background. 

The combinatorial background can be reduced by choosing a stau which gives the smaller 
value of the invariant mass M^r for every r candidate. Since the neutralinos produced by the 
decay of gluinos or squarks are likely to be highly boosted, the stau and r tend to be emitted in a 
similar direction, and therefore the invariant mass is likely to be much larger than the neutralino 
masses if we choose a wrong combination. To see how it works, we show the invariant mass 
distribution Mfr of the combination selected by the above strategy. We use the four-momentum 



of r extracted from the event generator (Fig. 11). Here, no information on the charge of r is 



used. As we see from the figure, the lowest mass combination shows peaks at the masses of 
Xi 2 3- The correct combinations are chosen with the probability of about 70%. 

In the actual experiment, however, the four-momentum of r-lepton is not available. Instead, 
the four-momentum of the r-jet from the hadronic r decays is available. If we use the four- 



momentum of the r-jet in the Mfj- analysis, the peaks shown in Fig. 11 are smeared by the 
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Figure 11: The f — r invariant mass distribution. The combination with the lowest invariant 
mass is chosen. The four-momentum of r-lepton extracted from the event generator is used. 



effect of the missing energy in the r decays. The left panel of Fig. 12 shows the distribution of 
the energy fraction of r-jet, Er-jet/Er, in the neutralino decays. We plotted histograms for each 
neutralino, Xi; X2 Xs- We rescaled the histograms so that the number of events are the same 
for each neutralino. Energies are measured in the laboratory frame. With the HERWIG event 
generator and the TAUOLA package, effects of the polarization of r are taken into account. 
As we can see, there are sharp edges in the distribution at Er-jct/Er = 1. Especially, the 
edge is sharper for Xi compared to X2 ^^'^ Xs- This can be understood as an effect of the 
polarization of r.0 Since the stau is mostly right-handed, the chirality of r from the neutralino 
decay is right-handed (left-handed) if the neutralino is gaugino-like (Higgsino-like) . By the V — A 
current structure of the weak interaction, neutrinos tend to be emitted in the opposite (same) 
direction to the r direction if r is right-handed (left-handed), and that makes the edge sharper 
(broader) With this structure, we can expect that the Mf^- distribution reconstructed with 
T-jet four-momentum shows sharp edges at three neutralino masses although the Higgsino edges 
become slightly weaker. 

In the right panel of Fig. |l^ we plotted the same quantity, E-r-jet/Er, from all the neu- 
tralino decays. The overall shape, monotonically increasing function and has a sharp edge at 
Er-jet/Er = 1, Can be understood from the distribution of ET-jet/Er in the two-body decays of 



*We thank L. Dixon for pointing out the possibiUty of having polarization effects. 
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-E'r-jet/ Ej Er-jet/Er 



Figure 12: The distribution of the r-jet energy fraction Er-jct/Er in the hadronic decay modes 
of r in supersymmetric cascade decays. In the left panel, we show the energy fractions for r's 
which originate from three species of neutralinos, Xij X2 ^^'^ /tsi respectively. They are rescaled 
so that the number of events are the same for three neutralinos. In the right panel, we did not 
distinguish the origin of r. Shaded histograms are the distribution of the energy fraction for the 
two-body decays, r ttz^, t ^ pu and r aiu assuming stable mesons. Energy calibration of 
the T-jets is performed by AcerDET. 

r^7rj/(ll%), T^pv {26%), t ^ mi^ (18%). (64) 

The p and oi mesons subsequently decay into two pious and three pious, respectively. Here the 
percentages of each mode denote the branching ratios. The branching ratio of the leptonic modes 
are 35%, and the other 10% comes from more than five-body decay modes or the modes with 
K mesons. When we ignore the width of the mesons, the energy fraction E^cson/ Er from each 
decay distributes uniformly between {m^^^^^/m^, 1) in the relativistic limit of r {E^- S> m,-). In 
the figure, we show the distributions of E^cson/ Er as shaded histograms. The energy fraction 
in other hadronic modes tends to pile up near the edge because of the kinematics of the many 
body final state. The distribution of E^cson/ E-r well resembles the properties of the distribution 
of Er-jet I Er ■ The thresholds at each meson mass are smeared by the effects of their finite decay 
widths (see for e.g., [pO| ). 

It is important to notice that the distribution of E-j-jet/Er has a tail in the unphysical region, 
E-r-jet/Er > 1. Thcse entries come mainly from the calibration of the r-jet energy used in the 
detector simulation. Especially, we should note that the distribution is slightly biased toward 
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Figure 13: Left) The distribution of the lowest invariant mass combination of fi and r-jet. The 
shaded histogram shows the events with a mis- identified r-jet which is simulated by assuming 
a mis-tagging probability of a non-r-labelled jet to be 1%. The small allows and dashed lines 
denote the input values of three neutralino masses. Three curves are fitting functions of three 
endpoints which correspond to the endpoints of x? 2 a from left to right, respectively. The third 
endpoint is statistically not very significant. Right) The same as the left figure but we assumed 
the mis-tagging probability to be 5% per a non-r-labelled jet. The endpoints of X2 X3 
visible whereas the significance of Xi events are reduced due to the shape of the background 
events. The bin size is 10 GeV in the right figure. 



the ET-jet/Er > 1 region. The detailed shape of the distribution, of course, depends on the 
actual algorithm for the calibration. We performed a fitting of the distribution around the peak 
with a smeared jagged function f{x), 



f{x;xo,a,Ci,C2) 
9ix]Xo) 



-00 V 2 



x';xo) 



■ exp 



vro"^ 

Ci X, (0 < X < xq), 
C2X, (xo < x). 



x;^ 

'2ct2 



(65) 
(66) 



with four fitting parameters, the position of the edge xq, the smearing factor a, and two slopes 
Ci and C2. The fitting gives the position of the edge xq to be xq = 1.049 ± 0.003, about five 
percent larger than unity. Since we will identify the position of the edge in the M^r distribution 
as the neutralino mass, the bias ends up with systematic errors of the mass measurement toward 
larger values. Therefore, in the actual analysis of the LHC data, we need to understand the 
shift of the edge location caused by the calibration of the r-jets energy. 

Having understood the edge structure of the Er-jct/Er distribution, we try to reconstruct 



the neutralino masses. Fig. 12 shows the distribution of the smaller invariant mass out of two 
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TUedge [GeV] 


am [GeV] 


Ci 


C2 




<r [GeV] 




194 ±2 


3±2 


0.93 ±0.35 


0.37 ±0.35 


19 ±20 


187 




279 ±3 


13 ±4 


0.33 ±0.14 


-0.014 ±0.086 


20± 17 


276 


xl 


314 ± 1 


1 ± 1 


0.066 ± 0.027 


-0.018 ±0.018 


8.3 ±4.5 


307 



Table 2: The fitting parameters of the function f{x — nif-^ ; medge — "ifi , Cmj Ci, C2) in Eq. (jd 
plus a constant oq around each edge. The final column shows the actual masses of Xi 2 3- 



possible combinations of fi and the r-jet. We took the mis-tagging probability to be 1% (left) 
and 5% (right). As we expected, we can clearly see the edge structures in the left panel. We 
can determine the masses of Xi 2 3 from the location of the edges. The tail of the distribution 
for Mrr ^ 350 GeV stems from the mis-tagging of r-jets (shaded histograms). In addition to 
the background from mis-tagging of r-jets, the histogram includes events with a fake stau from 
muons in supersymmetric events. The standard model background is assumed to be negligible 
with the selection cut discussed before [S8|. 



In the left figure, we can see the structure that the shape in Fig. 11 is smeared according to 
the Er-jet/Er distribution in Fig. 12. The edge structure of the third neutralino X3 is not very 
clear with the bin size of 5 GeV. In order to extract the masses of Xi 2 31 '^^ have fitted the each 



endpoints with the smeared jagged function f{x — m^^; medge ■ 



nii 



, cr^, Ci, C2) in Eq. (||) plus 



a constant (oq) which represents the background events around the edges. The results of the 
fitting are given in Table. |2[ 

By taking into account the above physics as well as systematic uncertainties such as depen- 
dence on the calibration algorithm for the r-jet momentum, we conclude that the masses of first 
two (possibly three) neutralinos can be measured with an accuracy of, at least, about 5% level. 

If we use a loose strategy for the identification of r, the background will significantly affects 



the edge structure. In the right panel of Fig. 13, we have used the mis-tagging probability to 
be 5% per a non-r-labelled jet. The edge structure is not significant for x? whereas we can see 
the edges of X2 ^nd Xs- (We changed the bin size to 10 GeV in this figure.) This situation will 
improve when we use a looser cut on px of r-jets such as pT > 20 GeV. A similar accuracy for 
the neutralino mass measurement is possible even in that case. 



4.3 Parameter determination 



It is straightforward to determine the model parameters (/i, Mmess; ^) from the measurement of 
TTif^ , m^o , and m^o . Performing x^ analysis would either give best fit values of these parameters 
or exclude the model. 
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Figure 14: Left) The stau mass m^^ as a function of Mmess for four values of the /i-parameter. 
The overall scale is set for M = 900 GeV. The dashed horizontal line corresponds to fi mass, 
nif^ = 116 GeV, at the benchmark point. The thick vertical line denotes the value of Mmess 
determined by assuming 5% precisions of ^ and M. Right) The pseudo-scalar Higgs boson mass 
niA as a function of M^ess for four values of the /i-parameter. The overall scale is set for M = 
900 GeV. The dashed horizontal line corresponds to the prediction of niA for Mmess = lO^'^ GeV. 
The thick vertical line denotes the value of Mmess determined from the stau mass measurement 
(see the left panel) . The arrow on the niA axes denotes the error of the prediction including the 
error AM. 



In most cases, a simpler analysis than the global fit is possible. First, by assuming that the 
model is correct, we can find that one of the two neutralinos we measured in the previous section 
should be Higgsino-like since their masses deviate from a GUT relation between Mi and M2. 
Secondly, we can neglect the tan (5 dependence in the neutralino masses. With a large value of 
tan/3 (see Fig. |^), corrections are of 0(1/ tan/?). Thus, the neutralino masses depend merely 
on the messenger scale Mmess- The parameters /i and M can be determined at the level of 5% 
from the measurement of two leading neutralino masses. If we can also measure the mass of X3) 
we can check the consistency of the GUT relation between Mi and M2 from the mass splitting 
between X2 and which provides a non-trivial check of GUT theories. 

We can then determine Mmess from the measured stau mass. We demonstrate in the left 



panel of Fig. 14 the determination of the messenger scale Mmess- Once we know the value of ^ 
and M, ruf^ can be calculated as a function of Mmess- Since we measure nif^ at a few permille 
level, and /i and M at 5% level, we can read off the corresponding value of Mmess from the 
figure. We find that the exponent of Mmess is determined with an accuracy of ±0.2. Therefore, 
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at the benchmark point, the parameters are determined with the precision of, 



A/i~20GeV, (67) 
AM ~ 50 GeV, (68) 
AlogioMmess~0.2. (69) 

Finally, once all the parameters are determined, we can make a prediction of any other 
physical quantities. The simplest test is the peak location of the Meg distribution which depends 



on the squark masses |91, 83, p2[. As a more non-trivial example, we show a prediction of the 



mass of the pseudo-scalar Higgs boson, m^, in the right panel of Fig. 14, We can predict a value 
of niA with an uncertainty of 5% level: 

AruA ~ 40 GeV, (70) 

around ruA = 765 GeV. At the LHC, the heavy Higgs boson H^/A^ will be discovered up to 



i^H/A — 800 GeV for tan /3 ~ 40 assuming and integrated luminosity 30 fb ^ [93|. Therefore, we 
can perform a non-trivial check of the model from the tha measurement. 

5 Conclusions 

In studies of LHC signals of supersymmetric theories, setting the model parameters is the first 
non-trivial task which needs to be done. At the LHC experiments, the main production process 
of the supersymmetric particles is a pair production of colored objects, such as a pair production 
of gluinos and squarks. It is thus essential to know how these particles decay. Also, since there are 
many kinds of particles, it is often the case that a measurement of supersymmetric parameters 
suffers from background processes which also come from supersymmetric events. In order to 
estimate the amount of the background, we need to set all of the parameters in the Lagrangian. 

There are, on the other hand, over 100 parameters in the Lagrangian of the MSSM. It is 
practically impossible to study every point in the 100 dimensional parameter space. Therefore, 
parametrizations such as the mSUGRA model and the "gauge mediation" model (the gauginos 
and scalar masses from the formula of the gauge mediation |5[ ^] and fi and B parameters as 
free parameters) have been proposed and used as standard benchmark models. The number of 
the parameters is reduced to be a few in these models. 

As a purpose of the study of generic signatures of the supersymmetric models and for 
development of methods to extract physical quantities, these parametrizations have played a 



significant role in studies of collider physics |91, 83, 32, 66 1. However, it is dangerous to rely too 
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much on these parametrizations. Interesting parameter regions in the MSSM can be precluded by 
assumptions made without theoretical motivations. Remember that they are simply convenient 
parametrizations of the more than 100 unknown parameters. 

The sweet spot scenario we have presented provides an example of a simple paramctrization 
of the MSSM (by only three parameters), and it is theoretically supported. It is the first example 
of such a simple paramctrization which has in background a well-defined closed framework of 
the supersymmetry breaking and mediation with no phenomenological/cosmological problems. 

There are new features in the collider signatures. For example, a relatively light Higgsino is 
preferred and that makes the decays of gluino into third generation quarks to be the dominant 
channel. Many numbers of 6-jets will show up in each supersymmetric event. Also, the light 
stau is predicted if the Higgsinos are light. If it is the NLSP, two charged tracks left by escaping 
staus and r-jets from the x° — fr decays can be used to reconstruct the three parameters 
in the Lagrangian as we have demonstrated. The model confirmation/exclusion is possible by 
measuring any other quantities such as the mass of the pseudo-scalar Higgs boson. 
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